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Phase transition and topology in 4d simplicial gravity 
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We present data indicating that the recent evidence for the phase transition being of first order does not result 
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1. Introduction 

When Euclidean Quantum Gravity is dis- 
cretized following the dynamical triangulation 
recipe, the Einstein-Hilbert action takes the form 



S = K4N4 — K2-/V2 



(1) 



where iVj is the number of i-simplices and (k 2i K4) 
are coupling constants. For a given topology the 
grand canonical partition function Zqc is defined 
by 



(«2)«4) 



(2) 



where C(T) is a symmetry factor and the sum 
goes over simplicial complexes T. The canonical 
partition function Zc is defined by an analogous 
expression, where N4 of T is kept fixed. One can 
investigate the critical behaviour of the theory by 
studying 



F{K2,Ni) = logZ c ( K2 ,7V 4 ) 



(3) 



and its derivatives in the limit N4 — > 00. One 
finds that the model has two phases. The de- 
termination of the exact nature of the transition 
between these two phases is obviously of utmost 
importance. 

In Ref. 0] it was found, in spherical topology, 
that the transition is of first order, contrary to 
the common belief. However, one can wonder 
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whether this finding does not reflect a deficiency 
of the simulation set-up. As it is an important 
question, we have decided to check it in more de- 
tail. The result of this investigation is presented 
in sect. 2 of this communication. In sect. 3 we 
also present some results on the behaviour of the 
model in other topologies. 

2. The dependence of the transition on the 
limitations of the algorithm 

In 4d there does not exist a set of ergodic lo- 
cal moves preserving the value of N±. In order to 
study the canonical ensemble the standard proce- 
dure is to carry out a grand-canonical simulation 
with a modified action, for example 



Smod = ~K 2 N 2 + 



V A f 



(4) 



where V4 is the value of N4 to be studied. All 
five local moves ergodic in the grand-canonical en- 
semble are employed. The constant k\ is chosen 
so that N4 fluctuates around V4: AN4 = 1/^/7. 
Measurements are performed every n'th time N4 
crosses V4, with an appropriate choice of n. In 
Ref. § we have used 7 = 0.005 and 0.0005. The 
first-order nature of the transition was concluded 
from finite size scaling and from the appearance 
of a double-peak structure at k 2 = 1.258 for our 
largest value of N4 = 32000. With 7 = 0.005 
(0.0005) one has A7V 4 = 14 (45) ) a rather mod- 
est fluctuation compared to 32000. However, the 
constraint limiting volume fluctuations may not 
be innocuous. Any two lattice configurations at 




N± = Va can be connected by a sequence of local 
moves, provided N4 is allowed to deviate enough 
from V4. When this condition is satisfied, the 
path of the system, in its time evolution and in 
the space of configurations, recurrently crosses 
the hyperplane N4 = V4. This defines the density 
of states, corresponding to the canonical ensem- 
ble. It is not obvious that this density is unaf- 
fected when volume fluctuations are constrained. 
Therefore, we have performed another set of sim- 
ulations at K2 — 1.258 and N4 = 32000, allowing 
for fluctuations AN4 ~ 10 4 . One cannot just set 
7 = 10~ 8 : the finite size corrections to F(k2,N4) 
lead to metastable states at different values of N4. 
Thus, instead of (|J) we choose another action 

Smod = ~^N 2 + U(N A ) , (5) 

where the external potential U is chosen so as to 
favour large excursions from the reference volume 
V4. The optimal choice is U = —F(k2,N4,) > but 
the free energy is unknown. Hence, we measured 
the derivative dF/dN4 for a few volumes, inter- 
polated it linearly and integrated over N4. The 
result can differ within the errors bars from the 
true free energy and this can make the simulation 
unstable. To insure stability we added a constant 
attractive force S = 0.0001 to the derivative and 
reflecting walls at 16000 and 64000. The long ex- 
cursions are evident in the history of N4 as shown 
in Fig.l. We have performed 1.5 x 10 6 sweeps 



and gathered 3 x 10° measurements, which were 
taken every 100'th crossing of V4. In Fig. 2 we 
show the corresponding A^o-distribution together 
with that of Ref . |l[] . They are in fact very sim- 
ilar and, in particular, both are bimodal. The 
first three cumulants of Nq distribution are now 
0.1966(4), 0.318(26) and -3.0(30) to be compared 
respectively with the values 0.1971(2), 0.316(8) 
and —6.1(21) given in Ref. [jlj. They agree within 
the statistical errors. In Fig. 3 we show the his- 
tory of the average size of minimal neck baby uni- 
verses, large in the elongated phase and small in 
the crumpled one, and the order of the most sin- 
gular point, which has the correlated opposite be- 
haviour. The existence of two states with tunnel- 
ing between them is very apparent. We conclude 
that the two state signal, i.e. the first-order na- 
ture of the transition is a genuine feature of the 
model. 

3. Dependence on topology 

We have investigated the dependence of the free 
energy on topology. More precisely, we have mea- 
sured 

1^:=!^ (6) 

and 

k{n 2 ,N4)^^ (7) 
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MEASUREMENTS 



Figure 3. The time history of the maximal vertex 
order and of the average minbu size. 



In the limit N4 — > oo they become intensive quan- 
tities and k(K2,N4 : ) — > «4 cr . We have simulated 
the topologies S i ,S 3 x S 1 and (S* 1 ) 4 . The mini- 
mum configurations have (N2MIN, N^min) equal 
(20, 6), (110, 44) and (1472,704) respectively. In 
the elongated phase the most probable move, ie 
barycentric subdivision, leads to 



^2 = §A 4 + c 



— N 2 MIN — 2~N4 MIN 



(8) 



and c's are 5,0 and 288 for 5 4 ,5 3 x S 1 and (S 1 ) 4 
respectively. One expects 

r-l , 



cN; 



k = K icr + (7 - 3)^4 1 + 



(9) 



with 7 = i. In Fig. 4 we show r as a function of 
I/N4 for K2 = 2.0. In fact, as one can see from 
figure 4 one gets a very good fit by using c from 
the minimum configurations. 

In the crumpled phase one expects 
-5 
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Figure 4. The average action density r vs. system 
size for various topologies. 



(11) 



where the coeeficients fk , f r are functions of K2 
only. At K2 = 0.0 i.e. in the crumpled phase we 
have made very extensive runs. We see no depen- 
dence on topology there, and the best fit gives 
5 = 0.49(12). This is in contrast to some earlier 
results, and consistent with S = 1/2. The latter 
value can be obtained again from a purely kine- 
matical argument, using the fact that the most 
probable configurations have N4 ~ Nq . As a con- 
clusion, the termodynamic limit is independent of 
topology, and the finite size corrections can be un- 
derstood from purely kinamatic arguments. An 
extended version of our results on topology de- 
pendence will be published in Ref. M. 
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